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Abstract 

It is shown that atom interferometry ahows for the construction of MIGO, the Matter- wave 
Interferometric Gravitational-wave Observatory. MIGOs of the same sensitivity as LIGO or LISA 
are expected to be orders of magnitude smaller than either one. A design for MIGO using crystalline 
diffraction gratings is introduced, and its sensitivity is calculated. 
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I. INTRODUCTION 



Roughly speaking, since an atom with mass m "weighs" much more than a photon with 
frequency uj, all other things being equal, matter- wave interferometers are expected to be 
more sensitive to inertial effects than light- wave interferometers by a factor of mc^/ hu ~ 10^° 
if atoms are used Because of this factor, it is expected that observatories constructed 
from matter-wave interferometers will be much more sensitive to gravitational waves (GWs) 
than those constructed with light-based interferometers. This expected increase in sensitiv- 
ity makes use of the quantum phase of an atomic deBroglie matter wave X^b = 2TTh/mv 
where v << c is the speed of an nonrelativistic object with mass m. When a GW passes 
through an atom interferometer, the speeds of the atoms in the interferometer — and hence 
their deBroglie wavelengths — will be slightly perturbed. By suitably designing the interfer- 
ometer, the shifts in the phase of the atoms caused by this perturbation in their deBroglie 
wavelength will result in a shift in the interference pattern, which can then be measured 
interferometrically. Like LIGO, this shift can be used to determine the properties of the 
GW. We shall argue in this paper that using matter-wave interferometry, a Matter-wave 
/nterferometric Gravitational-wave Observatory (MIGO) can be constructed. For the same 
sensitivity, MIGO is expected to not only be orders of magnitude smaller than laser-based ob- 
servatories such as LIGO (Laser interferometer Gravitational- wave Observatory) and LISA 
(Laser Interferometer iS'pace ylntenna), but is also expected to have a frequency response 
range much broader than current GW observatories. 

While classical-based systems such as LIGO and LISA place classical test masses (mir- 
rors) a certain distance way from the central beam splitter of the interferometer, and mea- 
sures slight shifts in their positions through light-based interferometry, quantum-based sys- 
tems such as MIGO throw out quantum test masses (atoms) with certain velocities, and 
measures the quantum phase shift due to slight shifts in their speeds through atom-based in- 
terferometry. As we will show below, given a MIGO with a characteristic length LmigOi for 
the same shot-noise-limited phase shifts as LIGO at 125 Hz (where LIGO is most sensitive) 



where 2B = 150 is the number of times light bounces back and forth in each of LIGO's 
arms, uj^ is the frequency of the laser used in LIGO, m is the mass of the atom used in 
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MIGO, / is the frequency of the GW, and Lmigo is the length of one of LIGO's arms. 
Thus, Lmigo ~ 1 m as compared to Lugo = 4 km if a cesium atom is used. Similarly, 
a comparison, for equal shot-noise- limited phase shift, of a space-based MIGO with LISA 
leads to a MIGO configuration that at least ten thousand times smaller than LISA. 

In this paper, we present one possible configuration for MIGO, and calculate the sensi- 
tivity of this configuration to GWs. We also outline the technologies that could potentially 
be used to construct it. Even though we restrict ourselves to those technologies that have 
already been demonstrated in various atom diffraction and interferometry experiments, our 
purpose in this paper is to explore the extreme limits of these technologies, and through 
these limits, establish the fundamental limitations to MIGO's capabilities. In this vein, we 
are particularly interested in crystalline diffraction gratings and mirrors. 

While the diffraction of helium atoms off of cleaved lithium fluoride crystals was demon- 
strated by Otto Stern over seven decades ago P|, and single-crystal Bonse-Hart neutron 
interferometers was used to measure gravitational effects three decades ago, atom inter- 
ferometers using the natural periodic structure of crystals as reflective diffraction gratings, 
i. e., crystalline gratings, to make beam splitters, have yet to be demonstrated. We will 
nevertheless consider only crystalline gratings in this paper. We shall show that MIGO's 
overall length depends crucially on the periodicity of the diffraction grating used in the beam 
splitter; the smaller the periodicity, the shorter the interferometer. While in principle any 
periodic structure can be used to construct diffraction gratings, the smallest of such possible 
structures is limited by the natural periodicity of crystals. With a periodicity on the order 
of Angstroms, diffraction gratings made from single crystals will have a periodicity orders of 
magnitude smaller than current nanofabricated gratings, whereas diffraction gratings based 
on light have a periodicity in the hundreds of nanometer range. MIGOs constructed with 
crystalline gratings will thus be the smallest matter-wave-based GW observatories that can 
conceivably be made. 

II. BRIEF REVIEW OF ATOM INTERFEROMETRY 

Matter-wave interferometry is based on the particle-wave duality principle of quantum 
mechanics. This principle states that every massive quantum object can have characteristics 
of both a particle and a wave. Which of these two characteristics it exhibits depends on 
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the properties the experimentahst wishes to measure or exploit. By exploiting the wave 
nature of particles, matter-wave interferometers have been constructed using "elementary" 
particles such as neutrons , as well as "composite" particles such as atoms 0, 0| , and Ceo 
fullerene molecules 0]. 

Like light-wave interferometry, matter-wave interferometers can be divided into three 
distinctive parts: the source emitting the interfering particle, the "atom optics" consisting 
of beam-splitters and mirrors, and the detector. Similar to the beams of light in light- 
wave interferometers, the beams of particles passing through the interferometer picks up 
a phase shift as well, but now A0 = (Sj/ — S^l) /h, and S'J/ and S'Jf are the action for a 
particle travelling along the two quasi- classical paths 71 and 72 along the two arms of the 
interferometer. Large throughputs of atoms will be needed for MIGO, as well as advanced 
"atom optics"; we focus therefore here on the atomic beam source, and the beam splitters 

Q 

and mirrors, and refer the reader to Scoles [7|] for an review of atomic-beam detectors. 



A. Supersonic Sources and 2D Optical Molasses CoUimation 

Supersonic sources were first developed in the 1960s by chemists to study of chemical 
n 

reactions ^2,], and they have a number of properties superior to the effusive sources (ovens) 

which are more often used by physicists. While the velocity distribution of atom beams 

from effusive sources is essentially a Maxwellian one that is fixed by the temperature of 

the oven, the velocity spreads of supersonic sources are much narrower. Fractional velocity 

spreads of 1% — 0.3% for helium — have been achieved from continuous supersonic sources, 

and even smaller fractional velocity spreads can be achieved using pulsed sources. Nearly 

n 

monoenergetic beams with very high intensities can thus be produced. Pritchard j^], for 
example, has produced 10^^ atoms/cm^ /s/sr sodium beams, 
very cold helium beams with 1.5 x 10^^ to 1.5 x 10^° atom/sr/s. 

Continuous Supersonic Sources: A typical continuous supersonic source functions as 
follows: A jet of gas from a high-pressure reservoir escapes supersonically in free expansion 
through a nozzle, consisting of a small orifice typically 10 to 100 microns in diameter, 
into a differentially pumped low-pressure chamber that has a larger orifice at its output 
end called the "skimmer". This skimmer has the appropriate geometry to skim away the 
hotter components of the outer perimeter of the rapidly expanding gas jet, thus leaving only 
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the intense, low-temperature central component of the atomic beam to enter into another 
differentially pumped chamber. Importantly, after the expansion the atoms in the beam no 
longer collide with one another. The beam is often collimated yet again using a second slit 
at the output end of yet another differentially pumped chamber before it enters the main 
vacuum chamber containing the atom-based device, such as an interferometer |^]. With 
successive stages of differential pumping by means of either diffusion or turbomolecular 
pumps, one can maintain an ultra-high vacuum in the main chamber that is often needed. 
In addition, using optical molasses techniques [ol, beam collimation can be achieved without 
throwing away most of the atomic beam, as is currently done when coUimating slits are used. 
The combination of 2D optical molasses collimation with effusive sources has been achieved 
by Kasevich |5]. However, its combination with supersonic sources, as we propose here, can 
yield the high-brightness atomic beams necessary for good signal-to-noise ratios in MIGO. 

Pulsed Supersonic Sources: Most pulsed supersonic sources function in much of the same 
way as continuous supersonic sources, but with the addition of a fast-acting valve to pulse the 
beam. Duty cycles of roughly 10~^ can be achieved resulting in very high fluences. Because 
of the short duty cycles, smaller pumps than those used in continuous supersonic sources 
can be used. Of particular interest is the pulse source developed by Powers et. al. [l^, which 



to generate the atomic beam. The 



uses laser ablation of metallic sources such as lithium 
high-pressure gas chamber is replaced by the laser-plus-metal assembly, and optical molasses 
can then be added to reduce the velocity spread of the pulsed supersonic beam as well. 



B. Atomic Beam Optics 

There are a number of approaches that can be taken to fabricate the "atom beam op- 



tics" — the beam splitters and mirrors — needed for the construction o 



MIGO. These optics 
and by also using 



have been constructed by using nanofabricated physical gratingp_P, ^ 
light through either standing waves or stimulated raman pulses [12]. We shall focus on the 
possibility of using cleaved surfaces of single crystals, which yield near-perfect crystalline 
surfaces, in the construction of optical elements for atomic beams. 

Atomic beam mirrors can be constructed either out of diffraction gratings, or from spec- 



ular reflection off of smooth surfaces. Diffraction orders of atomic 
surfaces have a history stretching back to Otto Stern's experiments 



Deams from crystalline 
and research in this 



area continues today [13| . Consequently, construction of zero-diffraction order, crystal-based 
mirrors is certainly possible. Of more interest for MIGO is the work of Hoist and Allison, 
however. They have demonstrated that the flexible, thin, hydrogen-passivated single-crystal 
silicon wafer surface Si(lll)-(1 x 1)H can be used as a specular mirror for focusing helium 
atomic beams The development of these mirrors can not only be used to construct 

collimating optics for the atomic beam, but they also demonstrated that planar, specular- 
reflection atomic mirrors constructed from single crystals are possible as well. We shall 
see that the sensitivity of MIGO is greater if planar, specular-reflection mirrors are used 
instead of diffraction gratings as reflection optics, and we will focus primarily on these types 
of mirrors here. 

Since the diffraction of atomic beams off of single crystals has been demonstrated, it 
is possible, in principle, to construct coherent beam splitters using single- crystal diffrac- 
tion gratings that work in reflection. Such diffraction gratings have a periodicity of a 
few Angstroms, and will be orders of magnitude smaller than that of current nanofabri- 
cated diffraction gratings. Alignment of interferometers constructed from crystal diffraction 
grating is much more challenging, however, and it is still an open question as to whether 
considerations arising from the Debye- Waller factor will allow for the construction of these 
interferometers. We address both issues below. 

Laser Interferometer Alignment: Since crystal-based atom interferometers must be 
aligned within an Angstrom or so, laser interferometry must be used to position and align the 
various optical elements of these interferometers. Although this is a difficult engineering task, 
there is no fundamental obstacle towards this goal. The sensitivity of the laser interferometry 
in distance measurements to make such measurements arises from the photon-shot noise limit 
5(f laser — ^ / {NphotonT'Y^'^ — lA/lyum ^ 10~^, where NphotonT is the number of photons in a 
given measurement time r interval. Hence we would need NphotonT — 10^ in order to perform 
these measurements. This is not difficult achieve with standard lasers. 

The argument in support of the feasibility of aligning crystal-based atom interferometers 
is further strengthened when one realizes what has already been achieved by LIGO. Unlike 
LIGO, where the mirrors of the interferometer are in free fall, i.e., placed as pendula on 
the end of kilometer-long interferometer arms, as we shall see below, the mirrors and beam 
splitters in MIGO can be fixed to an underlying, mechanically rigid frame. Furthermore, 



6 



the size of MIGO is on the order of meters, and not of kilometers as in the case of LIGO, or 
of milhons of kilometers as in the case of LISA. This should make the task of aligning MIGO 
orders of magnitude easier than that for LIGO. Indeed, the active-feedback- mechanism tech- 
nology developed for LIGO to stabilize the mirror-assembly-central-beam-splitter distance 
to within 10~^^ m could be transferred to MIGO as well jl5|J. 

Decoherence and the Debye- Waller factor: Quantum interference requires that there 
exists no possibility of "which-path" information for the helium atom inside the interferom- 
eter. This requires that diffraction of atomic beams off of beam splitters, and their reflection 
off of mirrors be elastic and coherent. (Scattering of atoms in the beam off of background 
gas in the ultra-high vacuum system is assumed to be negligible) . A measure of the inelastic 
versus elastic components of these processes is based on the Debye- Waller factor W in the 
intensity ratio I / Iq = exp{—2W), where / is the diffracted (or reflected) intensity and Iq is 
the incident intensity. The factor ly is a measure of the fluctuations of the atoms in the 
crystal that diffract the incident helium beam. For diffracted helium beams, W = B/a^, 



where a is the 
rule of thumb 



a 



attice constant of the crystal, and B ~ 0.5 at room temperature Q|. The 



is that W/12 < 0.1, for sharp, elastic diffraction patterns to be seen. The 
requirement for observing interference is more stringent, however; the probability of emitting 
even a single phonon during the diffraction process must be negligible, and W/12 < 0.01 is 



required 12|. For crystal diffraction grating and mirrors, W/12 ~ 10^'^ at room tempera- 
ture and decreases at lower temperatures. Thus it is highly probable that the zero-phonon 
process is the dominant one, and therefore quantum phase coherence is expected in the 
interferometer. The observation of high-visibility interference fringes for neutrons in the 



Bonse-Hart interferometer 



is evidence that this is indeed the case. 



III. LIGO, LISA, AND THE DETECTION OF GWS 

The great majority of the current experimental searches for GWs are based on laser 
interferometry. These detectors are scalable by design to such sizes that the detection of 
GWs becomes feasible. A number of research groups located throughout the world Q| are 
expecting to begin to collect data soon: GEO600, German-British collaboration; VIRGO, 
a French-Italian collaboration; TAMA300, a Japanese effort; and ACIGA, an Australian 
effort. The current US-based, international collaboration, is LIGO. In addition, a space- 

7 



based laser interferometer system LISA is currently in the initial planning stage. We will 
focus specifically on LIGO. 

LIGO is a set of three interferometers based at two locations separated by 3020 km: 
Hanford, Washington and Livingston, Louisiana [15]. All three instruments are based on 
Michelson interferometers with Fabry-Perot arms. The physical length of the each arm of 
the main LIGO interferometer is 4 km, and with the Fabry-Perot interferometers in each 
arm, the optical path of the arm is increased by the bounce factor, 2B = 150. At the end 
of each arm is a massive mirror suspended vertically as pendula within a vacuum chamber, 
and the location of the mirror assembly must be held in position within 10^^° to 10^^'^ m 
with respect to the center of the interferometer. 

In its current LIGO I configuration, LIGO is designed to detect GWs in the frequency 
range of 40 to 10^ Hz. At frequencies below 40 Hz seismic noise causes a rapid decrease in 
LIGO's sensitivity, and at frequencies above 125 Hz the shot-noise limit^^ begins to gradually 
limit LIGO's sensitivity as well. Between 40 and 125 Hz LIGO's sensitivity is limited by 
thermal noise in the suspension system of the end mirrors. Construction of LIGO I began 
in 1996, and the main interferometers were commissioned in 2001. The first science runs 
were started in June of 2002; data from these runs are currently being analyzed. Upgrading 
LIGO I to LIGO II, which is designed for GW astronomy, is expected to begin in 2006 and 
to be completed by 2007. 



IV. THE MATTER- WAVE /NTERFEROMETRIC GRAVITATIONAL- WAVE 
OBSERVATORY (MIGO) 

A. The Vertical MIGO Configuration 

Fig. ^is a diagram of the vertical MIGO configuration, one of the possible configurations 
for MIGO. We start here with a simplified description of the interferometer. Detailed analysis 
of the response of this configuration to the passage of GWs will be presented below. 

In the vertical configuration for MIGO shown in Fig. ^ alkali atoms are emitted from a 
pulsed supersonic source, and are then slowed down using Zeeman slower to velocities so 
slow that the acceleration due to gravity g dominates the trajectories of the atoms through 
the interferometer. Immediately after the supersonic source, the beam is then collimated 
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using 2D optical molasses to further narrow the transverse velocity spread of the atomic 
beam. After the beam is diffracted off the initial beam splitter atoms in the n = ±1 
orders, mirrors, which retroreflect the almost parabolic upwards trajectories, are placed at 
the top of these trajectories. (Atoms in the higher diffraction orders are not used.) The 
atoms are then reflected off of these mirrors and fall back downwards, almost retracing their 
almost-parabolic upwards trajectories. The left and right paths of the falling atoms are then 
recombined by the same beam splitter, and subsequently detected. The bilateral symmetry 
of this interferometer along the vertical line ensures that the gross motion of the atoms on 
the left side of the interferometer is reflected on the right side. The parabolic trajectories 
of the atoms is shifted slightly during the passage of a GW, and, as can be seen from the 
force lines in Fig. ^ will be slightly different for atoms travelling along the left hand path 
compared to the right hand path. It is this asymmetrical shift of the atomic trajectories by 
a X -polarized GW that leads to the overall phase shift. 



B. MIGO and LIGO: The Underlying Physics 

In this section, we focus on the physics that underlies both MIGO and LIGO, and compare 
their relative sensitivities to GWs from astrophysical sources. We shall emphasize only the 
physics here, and our arguments will be more physical than formal; an exact calculation of 
the MIGO phase shift can be found in the next section. As such, in this section we consider 
only the effects of a GW passing over freely-falling particles, and neglect all other external 
forces on the system. 



As is we 
system 
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1 known, when a GW passes through a system, its effects on all parts of the 
2ol | must be included. It is not possible to shield against gravitational forces as 
it is for the other forces of nature. Only differences in the motion between objects — such as 
the mirrors and the beam splitters in an interferometer — can be measured; their absolute 
motion cannot. As a consequence, the motion of an object relative to an observer is governed 
by the geodesic deviation equation j3] containing an effective tidal acceleration acting on 
the object arising from the local curvature of the spacetime. In the long wavelength limit, 
this equation reduces to 

d X \ d hij j . . 

where by construction the origin of coordinate system is placed at the observer at the beam 
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splitter. The metric deviations from flat spacetime hij{t), i. e., the amplitude of the GW, is 
a 3 X 3 tensor, and not a vector as it is for electromagnetic waves, although GWs have two 
polarizations — the + polarization and the x polarization — ^just like electromagnetic waves. 

The geodesic deviation equation lies at the heart of both LIGO and MIGO: LIGO on the 
classical level, and MIGO on the quantum level. For LIGO, one considers a test mass (a 
mirror) at a distance L' from the observer (the central beam splitter), which is taken to be 
the origin of the coordinate system. When a GW passes over the apparatus, it causes small 
shifts ^* in this distance so that since = + ^\ where |^*| << 

To lowest order, 

1 . /■* /■*' d'^h - 1 • 
^ -L^ / dt' / dt"^ = -UK,, 

where for t < 0, we assume hij = 0. Thus, in LIGO one puts a mass (a mirror) at t = a 
certain distance away from the observer, and observes fluctuations in its position due to 
the passage of a GW. These fluctuations are measured by light-based interferometry. The 
freely-falling end mirrors of LIGO are designed to be these classical test masses, and the 
phase shift due to changes in the length of its arms is Acp^jGO = '^t^ALligo/\i where 
is the wavelength of the laser, and ALugo is the optical path length difference of LIGO 
induced by the GW. For the two arms of the interferometer this difference for a + polarization 
GW with amplitude /i+ is ALugo = BLLiGoh+ where the factor of B accounts for multiple 
reflections of the light beam within each arm of LIGO within the Fabry-Perot interferometers 
Q- Thus, 

^(i>LIGO = -T-BLLiGoh+, (3) 

and depends only on the amplitude of the GW, and the length Lligo of LIGO's arms. 
Strictly speaking this expression for Acp^iGO is only applicable for GWs with frequencies 
~ 125 Hz; at 4 km, the arms of LIGO is roughly 1/6 the wavelength of a GW in the upper end 
of its frequency response spectrum, and causality has to be taken into account 21]. Although 
design characteristics of LIGO such as power recycling and shot-noise modifications by the 
Fabry Perot arms will change the sensitivity from this simple form, it is adequate for our 
purposes. Note that in LIGO the photons used in making the measurement are not the test 
particles being acted on by the GW; only the end mirrors of the interferometer are being 
acted on. 
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FIG. 1: Schematic of vertical MIGO configuration, with a x polarization GW incident normal to 
the plane of the interferometer. 

In contrast with LIGO which places an end mirror at a certain distance away and measures 
shifts in its position, in MIGO one throws out a test mass, i. e., an atom, at t = with a 
velocity Vi from the origin, and observes deviations in its velocity, and thus in its path in 
spacetime, due to the passage of a GW. To see this, we take Vi = Vi + Pi, where « \Vi\. 
Then from Eq. 

dt ^ 2 dt"^ 

Once again to lowest order, 

1 /■* d'^h- 1 
A ~ -V' / t'dt'^ ~ -VHh,, 
^ 2 Jq dt^ 2 ^ 

where hij is the rate of change of the strain field of the GW. The parallelism between the 
displacement ^* in LIGO, and the velocity shift /3j in MIGO, is readily apparent. Classically, 
the velocity shift (3i is extremely small, and virtually impossible to measure. Quantum 
mechanically, however, this velocity shift can be measured by means of quantum interference. 
A quantum particle such as an atom in a matter-wave interferometer has a deBroglie 
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wavelength XdB = 27Th/mv, where m is the mass of the atom and v is its speed. Changes 
in the atom's speed Av due to GWs resuh in changes to its local deBroglie wavelength, 
and hence its quantum phase. If T is the time it takes the atom to traverse MIGO, and if 
Lmigo is the characteristic size of MIGO, then \ Vi\ ~ Lmigo/T, and the shift in the atoms 
position ALmigo after it traverses the interferometer is ALmigo ~ LMiGoThy^/2 for a x 
polarization GW with amplitude hx- Like LIGO, this results in a phase shift, but now of 
the form 

47r TTl 

A(f)MiGO ~ T ALmigo = 2vr— L^^^q/Zix (4) 

where / is the frequency of the GW, and the deBroglie wavelength of the atom is now used 
instead of the wavelength of light. Thus, in contrast to LIGO, MIGO's phase shift A0a//go 
is proportional to the rate of change of the GW amplitude. Although simplistic, this rough 
derivation of the MIGO phase shift nevertheless elucidates the underlying physics of MIGO. 
The above expression for the MIGO phase shift is close to the exact equation we shall 
derive in the next section. Equation follows directly from Eqs. and Q by setting 
AcpMiGO = A(f)LiGO, and the amplitudes of the two polarizations of the GW equal. 

The relative sizes of MIGO compared with LIGO given in Eq. would seem to be 
counterintuitive. MIGO makes use of slowly-moving, nonrelativistic atoms to make its 
measurements, while LIGO would seem to make use of fast-moving photons. At first glance 
it would seem that MIGO should be less sensitive to GWs than LIGO by some power of 
Vatom/c. This, howcvcr, would be an erroneous argument. As outlined in the above, it is not 
the effect of the GW on the photons used in the laser interferometer that is being measured 
in LIGO; it is the effect of the GW on the test masses (the mirrors) which are at rest that 
is being measured. 

C. Phase Shift Calculation 

In this section we present an explicit calculation of the vertical MIGO phase shift. As 
shown in Fig. Q the origin of our coordinate system is place at the base of the beam splitter. 
In the absence of a GW the atoms travel along parabolic trajectories. If a GW with hij{t) 
is incident along the y-axis, these parabolic paths will be slightly perturbed, and we can 
calculate the resulting phase shifts for the atoms by taking x{t) = xo(t) + xi(t) where xo{t) 
is the path of the atoms in the absence of a GW, and xi{t) is the first-order perturbation 
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to this path. 

If the rth atom is sent through the beam sphtter at t = t^, the path of the atom in 
the n = +1 diffraction order in the absence of a GW is: Vq^ = v±, Vqz = v\\ — g(t — tr), 
Xq = v±{t—tr), and Zq = v\\ {t—t.f)—g{t—t^.Y/2 for tj. <t < t^+T /2, where T/2 is the time for 
the atom to travel from the beam sphtter to the mirror, and g is the local acceleration due to 
Earth's gravity. Since our focus is on GWs, we approximate g a.sa. constant, and neglect local 
curvature effects from stationary sources such as the Earth. For tj.+T /2 <t < t^+T. during 
which the atom travels from the mirror back to the beam splitter, its path is = —v±, 
vqz = v\\ —g{t — tr), xo = L±_ —v±(t — tr), and zq = v\\{t — tr) — g{t — trY 12, where L± = v±T 
is the width of the interferometer and L\\ = gL']_/8v'^ is its height. The velocities v± and 
f|l are the initial velocities of the atom as it leaves the diffraction-grating beam splitter, 
and f_L = 27ch/ma where a is the periodicity of the grating; v\\ = (ff — f^)^/^ where Vs is 
the velocity of the beam incident the beam splitter. A similar set of equations hold for the 
n = —1 order with v± — f±. 

The perturbed path xi satisfies the geodesic deviation equations of motion in Eq. ((21), but 
with X approximated as Xq on the right hand side. Like the unperturbed system, we solve 
the equations of motion for Xi by dividing the trajectory of the atom into two parts: first 
from the beam splitter to then mirrors, and then from the mirrors back to the beam splitter. 
The effect of the mirrors on the paths of the atoms are taken into account through the mirror 
boundary conditions Vizitr + T~ /2) = viz{t,r + T^ /2) and Vixit^ + T^ /2) = —vix(tr + T^/2); 
the mirrors are taken to be rigidly mounted on the frame of the interferometer. The — and + 
superscripts denote the atom approaching and leaving the mirror, respectively. This jump 
boundary condition on the velocity of the atoms along the x-axis is due to the instantaneous, 
impulsive force the mirror exerts on the atom when it hits the mirror. 

The action for the geodesic deviation equation can be written as 



Writing = Vq + vI, we find that for the n = +1 path the change in the atom's action S'^ 
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ion is 



due to the GW relative to the unperturbed action 

~ [vUU + T-) - Vi,{tr + T+)] + lL^v^K,{tr + T/2) 

m / -'2 



+L\\v^K,{tr + T/2) - - / dt {K^{t)vl + h,{t)vl{t)} 

( rtr+T/2 rtr+T 

-v±< / dt K^{t)voz{t) - / dt Kz{t)voz{t) 



tr+T/2 



1 



tr+T 



+ ^9 / dt [K,{t)xoit) + Kz{t)zo{t) - 2z,{t)] , 

after successive integration by parts, and using the geodesic deviation equation of motion. 
For the n = —1 path, we get S'_i{tr) from S'^i{tr) by taking v± —v± in the above. Then, 

-A(j)MiGO ~ Ls. [vt^itr + T-/2) + v^,j{tr + T-/2)] + 
m '- -' 

rtr+T/2 

2vMUU + T/2)-2vA / dt h,z{t)voz{t) 



dt Kz{t)voz{t) \ + 9 dth^zXoit) 

'tr+T/2 J Jtr 

-g dt [zt\t) - z^\t)\ , 

Jtr 

where the superscripts ±1 denote the trajectories for the ±1 orders, and we have used the 
fact that x^^it) = — Xq ^(t), while z^^it) = ZQ^(t). Then, 

ptr+T/2 

(vt^ + Vi,^)itr + T-/2) = L\\h^z{tr + T/2) + V\\h.,z{tr) '9 dt h,z{t), 

Jtr 

from Eq. (0), while for tr < t < + T/2, 

{zt^ - Z^^){t) = Xo{t) [h,z{t) + h^tr)] - 2V^ [ dt' h,z{t'), 

Jtr 

and for tr + T/2 <t <tr + T, 

{Zt' - Z^'){t) = L^Kzit) -V^[t- tr] {Kz{t) - Kz{tr)} 

-2v±Kz{tr + T/2)[t -tr- T/2] 

ft ftr+T/2 



/ rt Ptr + l/Z 

+2v^ / dt' Kz{t') - / dt' Kz{t') 

\ J tr+T/2 Jtr 



'tr+T/2 

Integrating, and then taking the Fourier transform, we get 

A^MiGoif) = -27r^LxL||2//ix(/)e-*'^^'^T,(/T/4), 
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for a GW with amplitude hx{f) and frequency /. Only the x polarization causes a net 
phase shift in the interferometer; the + polarization does not, due to the bilateral symmetry 
of the interferometer. The resonance factor is 

F,(/T/4) = 1 - {sinc(7r/T/2)}2 + — ^ {1 - sinc(7r/T)} 

where sinc(x) = sinx/x. Note that F^, — 1 as f'T oo, and thus MIGO is most sensitive to 
GW in the high frequency domain. This sensitivity is a direct result of the specular reflection 
mirrors, and their effects on the atoms are embodied by the jump boundary conditions that 
we used. If diffraction-grating mirrors were used instead, the high-frequency sensitivity of 
MIGO would be drastically reduced. 

While A0j\//Go oc L±L\\, because the acceleration due to gravity, L\\ = gL\/Sv\, where 
fx, the horizontal velocity of atom after the beam splitter, is inversely proportional to the 
periodicity of the diffraction grating used. In particular, we see that L\\ oc o?\ the smaller 
the periodicity of the diffraction grating, the shorter the interferometer; hence the use of 
crystalline diffraction gratings. We thus see AcpMiGO for the vertical MIGO configuration 
is proportional to L\, and not to L± as it is for LIGO. Indeed, taking the absolute value of 
A'pMiGO) the shot-noise-limited sensitivity for the vertical MIGO is 



27rmLxL||/|F,(/r/4)|iVi/2 

where N is the number of atoms passing though the interferometer per second. In compar- 
ison, h{f)shot'^ = {fiuJr^/2Io'qcY^'^ (27r//co'^), where /„ is the power of the laser, and t] is the 



photodetector efficiency (Eq. 123aof Thorne 2j|). Thus, while h{f)^l^^ decreases at higher 



frequencies, ^(/)fiof'^ increases. This complementarity between the two sensitivities is due 
to fundamental differences in the signals measured by LIGO and MIGO: LIGO measures 
position, while MIGO measures velocity. 

In the above analysis the mirrors of the interferometer is assumed to be rigidly attached 
an infinitely-rigid interferometer that does not move when the GW passes through. This 
is an approximation. A more accurate model is to attach each mirror to a spring with a 
resonance frequency /o, which depends on the size of the interferometer as well as its material 
properties, and a quality factor Q. When the GW passes through the interferometer it will 
induce a velocity shift in the mirror, which will in turn induce an additional shift in the 
atom's velocity. The net affect is to subtract the term —pl{P ~ fo + ^/o//Q)/2 from F^. 
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FIG. 2: The ability of MIGO to detect GWs from various classes of burst sources is plotted, and 
compared with that of LIGO I, LIGO II, and LISA. 

The small dips at the high-frequency end of the MIGO graphs shown in Fig. |21 are due to 
the shift in the mirror's position by the GW where a speed of sound of 6000 m/s was taken 
for the material used in constructing MIGO, and the mounting of the mirror to the frame 
has a. Q = 1. 



D. Potential Sensitivity 

The graphs for MIGO are plotted using hMicOi and = 10^^ atoms per second. An 
atomic beam of ^'^'^Cs atoms from a 2000 K pulsed supersonic source is slowed down using 
a Zeeman slower, and additional 2D optical cooling is used to narrow the velocity spread of 
the atoms perpendicular to the beam. A LiF crystal, which has a lattice constant of 2 A, is 
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used reflection diffraction grating for the beam splitter; the mirrors are constructed out of 
single LiF crystals as well. Two sizes of MIGOs are graphed. The flrst has a width L± = 35 
m and a height L\\ = 6.7 m; the atoms in the beam of this MIGO must be slowed by a factor 
of 69 by the Zeeman slower. The second has a width Lj_ = 240 m and a height = 316 m; 
the atoms in the beam must be slowed by a factor of 10 by the Zeeman slower in this case. 



Graphs of GW signals from burst sources in the flgure were replotted from Thome's 



22 



Fig. 9.4. The plots LIGO I, LIGO II and LISA have been updated from Thome's original 



221, Fritschel, and Cutler 



23[, respectively, using 



figure, and were generated from Thorne 
Eq. Ill of Thorne 

At 35 m X 6.7 m, the smaller of the two MIGO configurations graphed has the same 
sensitivity as LIGO II, in LIGO II's operating regime, while being a very small fraction of 
its size. Similarly, at 250 m x 316 km, the larger of the two MIGO configurations has the 
same sensitivity as LISA in LISA's operating regime, even though its length is only one ten- 
thousandths of LISA'S length of 5 million km. Importantly, in contrast to LISA, this MIGO 
is earth-based rather than space-based. Both MIGO configurations are sensitive to a wider 
range of GW frequencies than LIGO and LISA, and they would extend the reach of GW 
observatories into frequency ranges — in particular between 0.1 and 10 Hz — not currently 
accessible by either LIGO or LISA. Note also the sensitivity to high-frequency GWs of all 
three MIGO configurations graphed in Fig. 3, compared with the gradual loss of sensitivity 
of LIGO at frequencies above 125 Hz, and of LISA at frequencies above 2 x 10^^ Hz. This 
is due to differences in the signal being measured by matter-wave- and by light-wave-based 
interferometers: MIGO measures shifts in the atom's velocity, and thus the rate of change 
of the GW amplitude with time; the larger the velocity shift, the larger the phase shift. 
This rate of change increases at higher frequencies, and is the underlying reason why the 
smallest-amplitude GW that the three MIGO configurations can detect in Fig. 3 decreases 
as 1//^''^ at high frequencies. In contrast, LIGO (and LISA) measures the position of the 
end mirrors, and depends only on the amplitude of the GW. Indeed, because a GW is time 
varying, the smallest-amplitude GW that LIGO can detect must increase as f^^'^; it actually 
increases as /^^^ because of shot-noise limits to Acpuco- This difference between MIGO and 



LIGO can as be explicitly seen by comparing 



KffsiT mc'yN ) \2nLMiGO 
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in the shot-noise hmit where is the rate photons enter LIGO. 
V. CONCLUSION 

We have shown in this paper that in terms of sensitivity, size, capabihty, and flexibihty, 
the quantum methods embodied in MIGO have overwhelming advantages over the classical 
methods embodied in LIGO in detecting GWs. What before took many kilometers with 
LIGO could be done in a couple of meters, and what took millions of kilometers with LISA 
could be done with only hundreds of meters. Most importantly, many of the technologies 
needed in constructing a MIGO have already been separately demonstrated. Thus, quan- 



fundamental problems at the 



tum interference can become an important tool for study o 
intersection of quantum mechanics and general relativity [24 1 . 
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